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Abstract. We report on our on-going project to calculate proton decay matrix elements using 
domain-wall fermions on the lattice. By summarizing the history of the proton decay calculation 
on the lattice, we reveal the systematic errors of those calculations. Then we discuss our approach 
to tackle those uncertainties and show our prehminary results on the matrix elements. 
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Nucleon decay is one of the most important aspect that any (SUSY) GUT model has. 
>• I At low energy dimension-six operators are dominant contribution to the proton decay 

while higher dimensional operators are suppressed by the inverse power of the heavy 
mass (Mx). The dimension-six operators consist of three quark and one lepton fields. 
^ ! While the lepton part is treated trivially, the matrix element of the three-quark part 

O ; = £'j\u'^CPR/Ldj)PLu'' (1) 

between the initial proton and final pion (K or 7] meson) states receives a highly non- 
p il perturbative contribution from OCD, which we want to tackle in this study. The matrix 

^ ■ element has a tensor structure 



X 



n;p\i^/^.Jp;k) = Pi^[Wo-i4W,]up, (2) 



■ where q = k — p is the momentum transfer, Up is the proton spinor. The relevant form 

factor Wq is what we need since the ^ is practically zero by the on-shell condition of the 
outgoing lepton. 

The lattice gauge theory gives the first principle computational ground for the 
hadronic quantities like this matrix element of proton decay. In the first two calculations 
on the lattice |l2i|3|], the tree level chiral perturbation theory [4] was used to evaluate Wq 
from the low energy constant a and /3, 

aPtUp = {Ol^ijp), pPtUp = (Ol^tjp), (3) 

which are calculated on the lattice. This method is sometimes called the indirect method. 
Few years ago, JLQCD published on their large simulation of the nucleon decay matrix 
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element fl*], where they employed both direct and indirect methods. The direct method 
was first used by the authors of ref. |5]. However, the treatment of the form factors was 
improper, which led a large discrepancy in the results from direct and indirect methods. 
Once the direct method is treated properly, JLQCD found the discrepancy not so 
large, yet, 30 — 40% in most of the cases. Their results of the matrix elements are 3-5 
times larger than those from a model calculation commonly used, which pushes down 
the theoretical estimate of the life time of the proton, and makes much severe constraint 
on the GUT models. 

The existing calculations are all done with the Wilson fermion at a single lattice 
spacing. The Wilson fermion has an 0{a) discretization error, where a is the lattice 
spacing. There are two sources of error propagating to the matrix elements. One is the 
measurement of the matrix element in lattice unit. The other is the estimate of the lattice 
scale a^. Even for the state of the art calculation by JLQCD, the systematic error of 
the is as much as 30% ^. Of course there should be a scaling violation for Wq, 
a and /3, too, which could diminish the overall violation by compensating that from 
the scale. But it is unknown until it is studied. Also the Wilson fermion breaks chiral 
symmetry explicitly. Thus, the applicability of the chiral perturbation theory at a finite 
lattice spacing is not guaranteed. One has to take the continuum limit of quantities of 
interest. 

The second problem is that up to now the operator renormalization has been done by 
one-loop (tadpole-improved) perturbation theory. This should be improved by employ- 
ing a non-perturbative technique ^. 

Finally the calculations are all done in the quenched approximation, where all quark 
loop effects are neglected. This approximation is commonly used in the lattice calcula- 
tion as the unquenched simulation is much more expensive. One has to check how large 
is the effect of quenching by doing the unquenched simulation. 

We use the domain-wall fermions in our simulation. This fermion dis- 

cretization has almost exact chiral symmetry and exact flavor symmetry. Hence there 
practically is no mixing of the operators with different chiral structure, making the data 
cleaner. Moreover, there is no 0{a) discretization error. This second point has been 
demonstrated in the simulation results for the hadron spectrum lO and the kaon B 



parameter fl4|,115|]. The chiral symmetry can be further improved dramatically by im- 
proving the gauge action lfl6L .13]. We use DBW2 gauge action which reduce the residual 
chiral symmetry breaking by factor 100 from that for the Wilson gauge action at a typical 
lattice spacing lll3ll . 



The dimensionless quantity, the product of the Sommer scale 1 6, 7] and the p mass romp is about 30% off 
from its continuum Umit fs^ at the simulation point of JLQCD. Note that the decay width is proportional 
to the square of Wq (dimension two), or a and )3 (dimension three). 
^ For the recent summary of the non-perturbative renormalization on the lattice, see ||^. 
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FIGURE 1. Ratio R{t) for the relevant form factor Wq for (^l.^lp) at mj = ;7i2 = 0.04. 



We use the 16^ x 32 lattice with a ^ ~ 1.3 GeV ^. The direct method uses the ratio of 
the three- and two-point functions 



{Uh)jUt)){Jp{t)Jp{to)) 



•p-> 



(4) 



where the proton and pion interpolating fields are located at ?o = 6 and t\ = 24 respec- 
tively. Momentum ±p with pa = (l,0,0);r/8 or (1, l,0);r/8 is injected to the pion and 
the operator in the three point function, as well as in the pion two point function in the 
denominator. y/Z^ and y^Z^ are overlap of J^,; and Jp to the corresponding pion and 
proton states, which is estimated from the fit of two point functions. 

Figure [T] shows ratio at a parameter point with the particular projection and subtraction 
to get Wo, which is taken from the fit to the plateau. In addition to the data shown in the 
ref. Ill 7.1 . we have further performed the calculation for the non-degenerate quark mass 
mi, m2 in the final pseudoscalar state, where the initial proton state is made up of quarks 
with mi. Then we get Wq as a function of mi, m2, and q^. The chiral perturbation ^ 
helps to fit Wq to get to the physical point. The results for various decay amplitudes are 
shown in Fig. |2l We are assuming the SU(2) symmetry for the u and d quarks. There 
are other possible matrix elements, but they can be calculated with the matrix elements 
in the figure when the SU(2) symmetry is intact. We are yet to have the renormalization 
factor for the operators by a non-perturbative renormalization. Preliminary value using 
the perturbative estimate of the renormalization factor [18] is listed in ref. [IJ^. 

The results of the indirect method are also shown in Fig. |2| The direct and indirect 
calculations give consistent results within the error, in contrast to the result of JLQCD. 
However, this could be caused by larger statistical error in our calculation. We need to 
have more statistics to judge it. Nevertheless, the relative size of the matrix element in 
our calculation for each decay mode is similar to that obtained by JLQCD. 



The more precise description of our simulation is given in ref. O. 
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<7I°|(ud)LUjp> 
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FIGURE 2. Summary of the relevant form factor of the nucleon decay with both direct and indirect 
methods in lattice unit. 



We have investigated the proton decay matrix elements at a lattice cut off of ~ 
1.3 GeV with the domain- wall fermion in the quenched approximation. The direct 
and indirect methods give consistent results within our precision. The non-perturbative 
renormalization program 1 1^ 2^ is underway to get the continuum matrix elements. 



Also we are performing the two flavor dynamical domain-wall fermion simulation, 
which will give us an idea of the size of the quenching error. 

We thank RIKEN, Brookhaven National Laboratory and the U.S. Department of 
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